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Abstract. The shock reflection problem is one of the most important problems in math- 
ematical fluid dynamics, since this problem not only arises in many important physical 
situations but also is fundamental for the mathematical theory of multidimensional con- 
servation laws that is still largely incomplete. However, most of the fundamental issues 
for shock reflection have not been understood, including the regularity and transition 
of the different patterns of shock reflection configurations. Therefore, it is important 
to establish the regularity of solutions to shock reflection in order to understand fully 
the phenomena of shock reflection. On the other hand, for a regular refiection configu- 
ration, the potential flow governs the exact behavior of the solution in C^'^ across the 
pseudo-sonic circle even starting from the full Euler flow, that is, both of the nonlinear 
systems are actually the same in an physically signiflcant region near the pseudo-sonic 
circle; thus, it becomes essential to understand the optimal regularity of solutions for the 
potential flow across the pseudo-sonic circle (the transonic boundary from the elliptic to 
hyperbolic region) and at the point where the pseudo-sonic circle (the degenerate elliptic 
curve) meets the reflected shock (a free boundary connecting the elliptic to hyperbolic 
region). In this paper, we study the regularity of solutions to regular shock reflection 
for potential flow. In particular, we prove that the C^'^-regularity is optimal for the 
solution across the pseudo-sonic circle and at the point where the pseudo-sonic circle 
meets the reflected shock. We also obtain the C^'" regularity of the solution up to 
the pseudo-sonic circle in the pseudo-subsonic region. The problem involves two types 
of transonic flow: one is a continuous transition through the pseudo-sonic circle from 
the pseudo-supersonic region to the pseudo-subsonic region; the other a jump transition 
through the transonic shock as a free boundary from another pseudo-supersonic region 
to the pseudo-subsonic region. The techniques and ideas developed in this paper will be 
useful to other regularity problems for nonlinear degenerate equations involving similar 
difficulties. 



1. Introduction 

We are concerned with the regularity of global solutions to shock wave reflection by 
wedges. The shock reflection problem is one of the most important problems in mathematical 
fluid dynamics, which not only arises in many important physical situations but also is 
fundamental for the mathematical theory of multidimensional conservation laws that is still 
largely incomplete; its solutions are building blocks and asymptotic attractors of general 
solutions to the multidimensional Euler equations for compressible fluids (cf. Courant- 
Friedrichs [M], von Neumann [36j, Glimm-Majda [21], and Morawetz [33]; also see [2l l8l \20\ 
[221 EH [an [35]). 

In Chen-Feldman [TUj, the first global existence theory of shock reflection configurations 
for potential flow has been established when the wedge angle 6^ is large, which converge 
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to the unique solution of the normal reflection when 9^ tends to tt/2. However, most of 
the fundamental issues for shock reflection by wedges have not been understood, including 
the regularity and transition of the different patterns of shock reflection configurations. 
Therefore, it is important to establish the regularity of solutions to shock reflection in order 
to understand fully the phenomena of shock reflection, including the case of potential flow 
which is widely used in aerodynamics (cf. [51 [T31 [HJ [Ml 133] ) ■ On the other hand, for the 
regular reflection configuration as in Fig. 1, the potential flow governs the exact behavior of 
solutions in C^'^ across the pseudo-sonic {sonic, for short below) circle P1P4 even starting 
from the full Eulcr flow, that is, both of the nonlinear systems are actually the same in a 
physically significant region near the sonic circle; thus, it becomes essential to understand 
the optimal regularity of solutions for the potential flow across the sonic circle Pi Pi and at 
the point Pi where the sonic circle meets the reflected shock. 




Figure 1. Regular Reflection Configuration 



In this paper, we develop a mathematical approach in Sections 2-4 to establish the 
regularity of solutions to regular shock reflection with the configuration as in Fig. 1 for 
potential flow. In particular, we prove that the C^'^'^-regularity is optimal for the solution 
across the open part P1P4 of the sonic circle (the degenerate elliptic curve) and at the point 
Pi where the sonic circle meets the reflected shock (as a free boundary). The problem 
involves two types of transonic flow: one is a continuous transition through the sonic circle 
P1P4 from the pseudo-supersonic (supersonic, for short below) region (2) to the pseudo- 
subsonic (subsonic, for short below) region Q; the other is a jump transition through the 
transonic shock as a free boundary from the supersonic region (1) to the subsonic region 
ri. To achieve the optimal regularity, one of the main difficulties is that the part P1P4 of 
the sonic circle is the transonic boundary separating the elliptic region from the hyperbolic 
region. Near P1P4, the solution is governed by a nonlinear equation, whose main part has 
the form: 

{2x ~ atlj^)tlj^^ + bipyy ~ tp^ ^ 0, (1.1) 

where a,b > are constants, and which is elliptic in {x > 0} where ip > 0, with elliptic 
degeneracy at {x = 0} where ip = 0. We analyze the features of equations modeled by 
(jl.ip and prove the C^'" regularity of solutions of shock reflection problem in the elliptic 
region up to the open part P1P4 of the sonic circle. As a corollary, we establish that the 
C'l.i-regularity is actually optimal across the transonic boundary P1P4 from the elliptic to 
hyperbolic region. Since the reflected shock P1P2 is regarded as a free boundary connecting 
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the hyperbolic region (1) with the elhptic region Q for the nonhnear second-order equation of 
mixed type, another difficulty for the optimal regularity of the solution is that the point Pi is 
exactly the one where the degenerate elliptic curve P1P4 meets a transonic free boundary for 
the nonlinear partial differential equation of second order. As far as we know, this is the first 
optimal regularity result for solutions to a free boundary problem of nonlinear degenerate 
elliptic equations at the point where an elliptic degenerate curve meets the free boundary. 
To achieve this, we construct two sequences of points such that the corresponding sequences 
of values of ip^x have different limits at Pi; this is done by employing the C^'" regularity 
of the solution up to P1P4 excluding the point Pi, and by studying detailed features of the 
free boundary conditions on the free boundary P1P2, i.e., the Rankine-Hugoniot conditions. 

We note that the global theory of existence and regularity of regular reflection config- 
urations for the polytropic case 7 > 1, established in [10] and Sections 2-4, applies to 
the isothermal case 7 = 1 as well. The techniques and ideas developed in this paper will 
be useful to other regularity problems for nonlinear degenerate equations involving similar 
difficulties. 

The regularity for certain degenerate elliptic and parabolic equations has been studied 
(cf . m Uni [301 [33 [35] and the references cited therein) . The main feature that distinguishes 
equation (jl.ip from the equations in Daskalopoulos-Hamilton [TB] and Lin- Wang [3D] is the 
crucial role of the nonlinear term —aipx'4'xx- Indeed, if a = 0, then (jl.ip becomes a linear 
equation 

2xipxX +blpyy - TpX ^ 0. (1.2) 

Then ^po(x,y) := cx^^"^ is a solution of (|1.2|) . and V'o with c > also satisfies the conditions: 

V' > in {a: > 0}, (1.3) 

= on {a; = 0}, (1.4) 

dy^p = on{y = ±1}. (1.5) 

Let V be a solution of ([Il2]) in Qi := (0,1) x (-1,1) satisfying ([OJ-IIISl). Then the 
comparison principle implies that ip > ^0 in Qi for sufficiently small c > 0. It follows that 
the solutions of (|1.2p satisfying (|1.3p - (|1.5p are not even C^'^ up to {x = 0}. On the other 
hand, for the nonlinear equation (|l.ip with a > 0, the function tp{x, y) = -^x"^ is a smooth 
solution of (|l.ip up to {x — 0} satisfying p.3p - (|1.5p . More general C"^'^ solutions of (|l.ip 
satisfying (|1.3p - (|1.5p and the condition 

2-/3 

- Mx <il^x< -X in X > 0, (1.6) 

a 

with M > and /3 e (0, 1), which implies the ellipticity of (jl.ip . can be constructed by using 
the methods of [10 , and their C^'" regularity up to {x — 0} follows from this paper. Another 
feature of the present case is that, for solutions of (jl.ip satisfying (|1.3p - (|1.4p and (|1.6p . we 
compute explicitly D^ip on {x — 0} to find that tpxxi^, u) — ^ ^ind ipxyiO, y) = 4'yy{^: y) = 
for all such solutions. Thus, all the solutions are separated in C^'" from the solution = 0, 
although it is easy to construct a sequence of solutions of (|l.ip satisfying (|1.3p - (|1.4p and 
(|1.6p which converges to -0 = in C^'". This shows that the C^'^-regularity of pTT]) 
with conditions l|1.3p - (|1.4p and (|1.6p is a truly nonlinear phenomenon of degenerate elliptic 
equations. 

Some efforts have been also made mathematically for the shock refiection problem via 
simplified models, including the unsteady transonic small-disturbance (UTSD) equation (cf. 
Keller-Blank [26,, Hunter-Keller [25], Hunter [24], Morawetz [33]) and the pressure gradient 
equation or the nonlinear wave system (cf. Zheng [39], Canic-Keyfitz-Kim [6]). On the other 
hand, in order to understand the existence and regularity of solutions near the important 
physical points and for the reflection problem, some asymptotic methods have been also 
developed (cf. Lighthill [29], Keller-Blank [26], Hunter-Keller [25], Harabetian [23], and 
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Morawetz Also see Chen [Tl| for a linear approximation of shock reflection when 

the wedge angle is close to ^ and Serre [M] for an apriori analysis of solutions of shock 
reflection and related discussions in the context of the Euler equations for isentropic and 
adiabatic fluids. We remark that our regularity results for potential flow near the sonic 
circle confirm rigorously the asymptotic scalings used by Hunter-Keller [25j . Harabetian 
[23] . and Morawetz [33]. Indeed, the C^'" regularity up to the sonic circle away from Pi 
and its proof based on the comparison with an ordinary differential equation in the radial 
direction conflrm their asymptotic scaling via the differential equation in that region. The 
optimal C^'^ regularity at Pi shows that the asymptotic scaling does not work there, i.e., 
the angular derivatives become large, as stated in [33]. 

The organization of this paper is the following. In Section 2, we describe the shock re- 
flection problem by a wedge and its solution with regular reflection configuration when the 
wedge angle is suitably large. In Section 3, we establish a regularity theory for solutions near 
the degenerate boundary with Dirichlet data for a class of nonlinear degenerate elliptic equa- 
tions, in order to study the regularity of solutions to the regular reflection problem. Then 
we employ the regularity theory developed in Section 3 to establish the optimal regularity 
of solutions for 7 > 1 across the sonic circle P1P4 and at the point Pi where the sonic circle 
P1P4 meets the reflected shock P1P2 in Section 4. We also established the C^'" -regularity 
of solutions in the subsonic region up to the sonic circle P1P4. We further observe that the 
existence and regularity results for regular reflection configurations for the polytropic case 
7 > 1 apply to the isothermal case 7=1. 

We remark in passing that there may exist a global regular reflection configuration when 
state (2) is pseudo-subsonic, which is in a very narrow regime (see [T31[33]). In this case, 
the regularity of the solution behind the reflected shock is direct, and the main difficulty of 
elliptic degeneracy does not occur. Therefore, in this paper, we focus on the difficult case 
for the regularity problem when state (2) is pseudo-supersonic, which will be simply called 
a regular shock reflection conflguration, throughout this paper. 

2. Shock Reflection Problem and Regular Reflection Configurations 

In this section, we describe the shock reflection problem by a wedge and its solution with 
regular reflection conflguration when the wedge angle is suitably large. 

The Euler equations for potential flow consist of the conservation law of mass and the 
Bernoulli law for the density p and the velocity potential $: 

9tp-|-divx(pVx$) = 0, (2.1) 

5t$-l-i|Vx$P+i(p) = So, (2.2) 

where Bq is the Bernoulli constant determined by the incoming flow and/or boundary con- 
ditions, and 

i'{p)^p'{p)/p^c\p)/p 
with c{p) being the sound speed. For polytropic gas, by scaling, 

p(p)-pV7, c2(p) = p^-i, ^{p)^P^^-l, 7>1. (2.3) 

7-1 

2.1. Shock Reflection Problem. When a plane shock in the (x, t)-coordinates, x = 
{xi,X2) £ K^, with left-state (p, Vx*!") — (pi,ui,0) and right-state (po, 0, 0), Ui >0,pq < pi, 
hits a symmetric wedge 

W :— {(xi, X2) : |a;2| < a^i ta,n9w,xi > 0} 
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head on, it experiences a reflection-diffraction process. Then the Bernoulh law (|2.2p becomes 



This reflection problem can be formulated as the following mathematical problem. 



(2.4) 



Problem 1 (Initial-Boundary Value Problem). Seek a solution of the system of equations 
p.ip and ()2.4I) . the initial condition at t = 0: 



(P,*)|t=0 



(po,0) /or I X2 1 > tan 6*^, > 0, 



{pi,uixi) for xi < 0, 
and the slip boundary condition along the wedge boundary dW : 

V$ • v\aw = 0, 
where v is the exterior unit normal to dW (see Fig. 2). 



(2.5) 
(2.6) 




Figure 2 . Initial-Boundary Value Problem 



Notice that the initial-boundary value problem (|2.ip and (|2.4p - (|2.6p is invariant under 
the self-similar scaling: 

$ 

(x, f) (ax, at) , (p, <J>) ^ (p, — ) for a ^ 0. 

a 

Thus, we seek self-similar solutions with the form 

p(x,<) -p(e,r;), $(x,i) -tV(C,?y) for = 

Then the pseudo-potential function = "0 ^ 5 (C^ + '7^) is governed by the following potential 
flow equation of second order: 

div {p{\D^\'', ^)Dip) + 2p{\D^\'',ip) = 

with 



(2.7) 



1, 



p{\D^\\ip) = {pV - {l~l){^ + ^\Dip\^))—\ 



(2.8) 



where the divergence div and gradient D are with respect to the self-similar variables (^, 77). 
Then we have 



= c\\Dv\\ ^, pr') = pV - (7 - 1)(2 I^^P + V)- 



(2.9) 
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Equation (|2.7p is a nonlinear equation of mixed elliptic-hyperbolic type. It is elliptic if 
and only if 

\D^\<c{\D^\\^,pT'), (2.10) 

which is equivalent to 



\D^\ < c, [if, po, 7) y {pV - (7 - 1)^) ■ (2.11) 

Shocks are discontinuities in the pseudo-velocity Dtp. That is, if fl'^ and :— Q \ 0+ are 
two nonempty open subsets of C and S :— 9f2+ n f2 is a C^-curve where D(p has a 
jump, then if € wl^^in) n C^Q"^ U 5) n C2(f]±) is a global weak solution of ()2.7|) in H. if 
and only if ip is in W^^^ (fl) and satisfies equation p.7p in il* and the Rankine-Hugoniot 
condition on 5: 

[p(|i^^P,^)i?V'-^]s = 0. (2.12) 

The plane incident shock solution in the (x, i)-coordinates with states {p, Vx^*) = {po, 0, 0) 
and (pi,ui,0) corresponds to a continuous weak solution ip of (|2.7p in the self-similar coor- 
dinates (^, 77) with the following form: 

M^,v) = -lie+v') for O^o, (2.13) 



2 

Vi{tri)=-\ie+v') + M^-^o) for C<^o, (2.14) 



respectively, where 



, _ , . / 2(P? ' - Po ') _ Piui , . 

V (7- 1)(PI -Po) pi - po 
is the location of the incident shock, uniquely determined by {po, pi, 7) through (|2.12p . that 
is, Po = (Coi Co tari^^) in Fig. 1. Since the problem is symmetric with respect to the axis 
77 = 0, it suffices to consider the problem in the half-plane 77 > outside the half-wedge 

A := < 0,7/ > 0}U {77 > Ctan6'„,, ^ > 0}. 

Then the initial-boundary value problem (|2.ip and (|2.4p - (|2.6p in the (x, i)-coordinates can 
be formulated as the following boundary value problem in the self-similar coordinates (C, 77). 



Problem 2. (Boundary Value Problem) (see Fig. I). Seek a solution ip of equation (j2.7[) 
in the self-similar domain A with the slip boundary condition on the wedge boundary dA: 

D^-i^\oA^O (2.16) 

and the asymptotic boundary condition at infinity: 

jvo /or C > Co,?7 > Ctan6l„, e2, 2 /^1-7^ 

ip^ip:^i^ whenC''+V (2-17) 

[fi for 4 < 4o,?7 > 0, 

where (|2.17p holds in the sense that lim \\^ ^'ip\\c(A\Ba{o)) — 0- 

2.2. Existence of Regular Reflection Configurations. Since pi does not satisfy the 
slip boundary condition (|2.16p . the solution must differ from (pi in < ^o} H A and thus a 
shock diffraction by the wedge occurs. In TOJ, the existence of global solution ip to Problem 
2 has been established when the wedge angle 9w is large, and the corresponding structure of 
solution is as follows (see Fig. [1]): The vertical line is the incident shock S = = Co} that 
hits the wedge at the point Po = (Co, Co tan 0^,), and state (0) and state (1) ahead of and 
behind S are given by (po and ip>i defined in (|2.13p and (|2.14p , respectively. The solutions ip 
and (fii differ within {C < Co} only in the domain P0P1P2P3 because of shock diffraction by 
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the wedge vertex, where the curve PqPiP2 is the reflected shock with the straight segment 
PoPi- State (2) behind PqPi is of the form: 

'f2{£.,v) = + V^)+M£.~£.o) + (?y-Cotan6'„)u2tan6i^, (2.18) 

which satisfies 

Dip -1^ = on aAn > 0}; 

the constant velocity U2 and the angle between PqPi and the ^-axis are determined by 
{&w,Po, Pi, 7) from the two algebraic equations expressing (|2.12p and the continuous match- 
ing of ipi and ip2 across PqPi- Moreover, the constant density p2 of state (2) satisfies p2 > pi, 
and state (2) is supersonic at the point Pq. The solution ip is subsonic within the sonic circle 
for state (2) with center (^2, ""2 tan6'n,) and radius C2 = p^2^^^^^ > (the sonic speed of 
state (2)), and ip is supersonic outside this circle containing the arc P1P4 in Fig. [1] so that 
P2 is the unique solution in the domain PqPiP^, as argued in [51 134j. Then p differs from 
P2 in the domain = P1P2P3P4, where the equation is elliptic. 

Introduce the polar coordinates (r,6) with respect to the center (^2, 1*2 tan 0^) of the 
sonic circle of state (2), that is, 

^ — 1*2 = r cos 6*, 77 — U2 tan = r sin 6*. (2-19) 

Then, for e S (0, C2), we denote hy il^ :~ D {{r, 0) : < C2 — r < e} the e-neighborhood 
of the sonic circle P1P4 within fi. In fl^, we introduce the coordinates: 

x^C2-r, y = 0- 0.U,. (2.20) 

This imphes that fie C {0 < x < e, y > 0} and P1P4 C {x — y > 0}. Also we introduce 
the following notation for various parts of dQ: 

dnndBc^{{u2,U2ta.n0^)) ee P1P4; 
P1P2; 

dnndA = PsPi. 

Then the global theory established in lOj indicates that there exist 0c = 0c{po, Pi,j) G (0, ^) 
and a = a{po,pi,^) G (0, 1) such that, when 6^ e [6c, f ), there exists a global self-similar 
solution: 

$(x,i) =t^(^) + ^ for^eA, i>0, 

with 

p(x,i) = {pr' - (7 - im + ^ivx<i>n)^ 

of Problem 1 (equivalently. Problem 2) for shock reflection by the wedge, which satisfies 
that, for (^,77) = f , 

P e C"'i(A), 

P e C°°(f^)nCi'"(f2), 

f2 21^ 

{po for ^ > ^0 and 77 > ^tan^^, ^ ' ^ 

pi for C < Co and above the reflection shock P0P1P2, 
P2 in P0P1P4. 

Moreover, 

(i) equation (|2.7|) is elliptic in fJ; 

(ii) p > p2 in ft] 

(iii) the reflected shock P0P1P2 is at Pi and C°° except Pi, 



^ shock 
^ wedge 
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(iv) there exists Eq & (0, ^) such that ip e C^'^{n^g) n C^{n^g \ Tsonic)i moreover, in 
the coordinates (I2.20p . 

li^-¥'2||^Xfi^:- E sup(x^-+ = -2|a,^a^(^-^2)(x,y)|)<cx3; (2.22) 



0<fc+i<2 ^S^eo 

(v) there exists (5o > so that, in the coordinates (|2.20p . 

2 S 

Ida^i'p- 'P2)ix,y)\ < — -^x in Oeo; (2.23) 
7+1 

(vi) there exist uj > and a function y = f{x) such that, in the coordinates l|2.20p . 

^ea = {ix,y) : X e (0, eo), 0<y< /(x)}, 

(2.24) 

Tshock n {0 < X < eo} = {{x, y) ■■ X e (0, eo), y = fix)}, 

and 

ll/llci.M[o,eo]) < oo , ^>cc;>0 for 0<x<eo. (2.25) 

The existence of state (2) of the form (12.18^ with constant velocity (u2, U2 tan^^,), U2 > 0, 
and constant density p2 > pi, satisfying (|2.12|) and ipi = (p2 on PqPi, is shown in ^JOl Section 
3] for 9w S [6*0 f)- The existence of a solution ip of Problem 2, satisfying (|2.2ip and property 
(Irv)) follows from [TUl Main Theorem]. Property (|i]) follows from Lemma 5.2 and Proposition 
7.1 in [TU]. Property ([n]) follows from Proposition 7.1 and Section 9 in [TD], which assert 
that ip — if2 & ^, where the set K. defined by (5.15) in [TU], which implies property 
Property (|v|) follows from Propositions 8.1-8.2 and Section 9 in [10 . Property (|vi|) follows 
from (5.7) and (5.25)-(5.27) in [lOj and the fact that (p - ip2 e IC. 

These results have been extended in [11] to other wedge-angle cases. 



3. Regularity near the degenerate boundary for nonlinear degenerate 
elliptic equations of second order 

In order to study the regularity of solutions to the regular reflection problem, in this 
section we first study the regularity of solutions near a degenerate boundary for a class of 
nonlinear degenerate elliptic equations of second order. 

Wc adopt the following definitions for ellipticity and uniform ellipticity: Let C be 
open, u e C^(r2), and 

2 

Cu^ A^jix, Du)u^.j + B{x, Du), (3.1) 

where Aij (x,p) and B(x,p) are continuous on x R^. The operator L is elliptic with respect 
to u in if the coefficient matrix 

A{x,Du{x)) [A^j{x,Du{x))] 

is positive for every a; € fi. Furthermore, C is uniformly elliptic with respect to it in f2 if 

\I < A{x, Du{x)) < AI for every x e 

where A > A > are constants and / is the 2x2 identity matrix. 

The following standard comparison principle for the operator C follows from 19 , Theorem 
10.1]. 

Lemma 3.1. Let fl C be an open bounded set. Let u,v G C(0) n C^{^1) such that the 
operator C is elliptic in Q with respect to either u or v. Let Cu < Cv in Q and u > v on 
dfl. Then u > v in fl. 
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3.1. Nonlinear Degenerate Elliptic Equations and Regularity Theorem. We now 

study the regularity of positive solutions near the degenerate boundary with Dirichlet data 
for the class of nonlinear degenerate elliptic equations of the form: 

Ciil^ := {2x - aipx + Oi)-0a;x + 02i]-'xy + [b + 0:i)%l}yy - (1 + 0^)^-)^ + O^i/jy = in 

(3.2) 

^ > in Q+^, (3.3) 

i>^0 on 9Q+^ n {x = 0}, (3.4) 
where a,b > are constants and, for r, i? > 0, 

:= {(x, y) : xe (0, r), \y\ < R} C M^, (3.5) 
and the terms Oi{x,y),i — 1, . . . ,5 are continuously differentiable and 

\0^{x^^\OUx^^^ forfc = 2,...,5, (3.6) 

X X 

^^^'^""■■^^^ \DOkix,y)\<N for fc = 2,...,5, (3.7) 



X 

in {x > 0} for some constant N. 

Conditions (|3.6p - (|3.7p imply that the terms Oi,i — I,-- - ,5, are "small"; the precise 
meaning of which can be seen in Section [3] for the shock reflection problem below (also see 
the estimates in [lOj). Thus, the main terms of equation (j3.2p form the following equation: 

{2x - aipx)'>p3:x + b^pyy - ipx ^ inQ^j^. (3.8) 

Equation (|3.8p is elliptic with respect to ■(/' in {a; > 0} if ■(/'a: < In this paper, we 
consider the solutions that satisfy 

2 — 3 

- Mx <i)x< -X in Q+ (3.9) 

a 

for some constants M > and (3 E (0, 1). Then (|3.8p is uniformly elliptic in every subdomain 
{x > S} with S > 0. The same is true for equation (|3.2p in Q^j^ if r is sufficiently small. 

Remark 3.1. If f is sufficiently small, depending only on a, b, and N , then (|3.6p - (|3.7p and 
p.9p imply that equation (|3.2p is uniformly elliptic with respect to ip in ^ {x > 5} for 
any S e (0, |). We will always assume such a choice of f hereafter. 

Let ip £ C^^iQf r) a solution of p.2p satisfying (13. 9p . Remark [3.11 imphes that the 
interior regularity 

■4> £ C^-'"{QIb) for aU a £ (0, 1) (3.10) 
follows first from the linear elliptic theory in two-dimensions (cf. [19j Chapter 12]) to 
conclude the solution in C^'" which leads that the coefficient becomes C" and then from 
the Schauder theory to get the C^'" estimate (cf. [121 Chapter 6]), where we use the 
fact Oi £ C^{{x > 0}). Therefore, we focus on the regularity of t/j near the boundary 
{x ~ 0}r] dQfj^ where the elHpticity of p.2p degenerates. 

Theorem 3.1 (Regularity Theorem). Let a,b, AI, N, R > and (3 £ {Q,j) be constants. 
Let ip £ C{Qfj^) n C^iQfj^) satisfy ([33]) ([SJ]) . and equation ([321) in Q^r with 
Oi = Oi(x,y) satisfying Oi £ C^iOf^) and ((3J)) - ([SJ]) . Then 



tp £ C^^"(Qt^2,fl/2) V any a £ (0, 1), 
with ^ 

'4'xx[^,y) = i^xv{0,y) = il'yy{Q,y) = for all \y\< 
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To prove Theorem 13. 1[ it suffices to show that, for any given a e (0, 1), 

V' e C-^^^IqJ^) for some r e (0, f/2), (3.11) 

since ip belongs to C2'"(Q+2 ^.^^ C^ {x > r/2}) by ([3A0| . 
Note that, by jSSll-dS^l) and ([31]), it follows that 

0<^{x,y)<'^^x^ for all (a:,y)eQt^. (3.12) 

The essential part of the proof of Theorem 13.11 is to show that, if a solution -ip satisfies 
p.l2p . then, for any given a G (0, 1), there exists r e (0, f/2] such that 

|V'(x, v) - ^^'1 < Cx^^-^ for aU (x, y) G Q+^j^/g. (3.13) 

Notice that, although ijj^'^^ = is a solution of p.2p . it satisfies neither p.l3p nor the 
conclusion V'I'k (0, y) = ^ of Theorem 13.11 Thus it is necessary to improve first the lower 
bound of in p.l2p to separate our solution from the trivial solution ?/;'^"^ = 0. 

3.2. Quadratic Lower Bound of if). By Remark |3.1[ equation p.2p is uniformly elliptic 
with respect to ip inside Qf Thus, our idea is to construct a positive subsolution of p.2|) . 
which provides our desired lower bound of tp. 

Proposition 3.1. Let satisfy the assumptions in Theorem 13.11 Then there exist r G 
(0,f/2] and /i > 0, depending only on a,b, N, R,r, P, and infQ+ r]{x>f/2} such that 

i^{x,y) > tix^ on Qt^^^R,/iQ- 

Proof. In this proof, all the constants below depend only on the data, i.e., a, &, M, N, R, f, (3, 
and infQ+ n{2;>f/2} unless otherwise is stated. 
Fix yo with \yo\ < We now prove that 

'P{x,yf))> pix^ forxe(0,r). (3.14) 

We first note that, without loss of generality, we may assume that R = 2 and yo — 0. Oth- 
erwise, we set -ipix, y) := ip{x, yo + ^y) for all (x, y) G Q^fj- Then e CiQf^) ^ (^^(Q^f^) 
satisfies equation l|3.2p with (|3.6p and conditions (|3.3p - p.4p and (|3.9p in Qf2^ '^ith some 
modified constants a, 6, TV, and functions Oi, depending only on the corresponding quan- 
tities in the original equation and on R. Moreover, infQ+ r]{x>f/2} ^ ~ ^^^Q+ n{x>f/2} ^■ 

Then (|3.14p for tp follows from (|3.14p for -0 with yo — and R = 2. Thus we will keep the 
original notation with yo = and i? = 2. Then it suffices to prove 

V'(x,0) > fxx'^ for X e (0,r). (3.15) 

By Remark l3. II and the Harnack inequality, we conclude that, for any r £ (0,f/2), there 
exists a = a{r) > depending only on r and the data a, 6, N, R, f, and infQ+ r]{x>f/2} V'j 
such that 

■0 > cr on Q^3/2 n{x> r}. (3.16) 

Let r e (0,f/2), A: > 0, and 

o'(r) 

0</i<^ (3.17) 

to be chosen. Set 

w{x,y) -.^ fix^{l-y^)^kxy\ (3.18) 
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Then, using (|3.16p - (|3.17p . we obtain that, for aU x € (0, r) and \y\ < 1, 

'w(0,y) = 0< V(0,2;), 

^w{x,±l) = -kx < < V(a;,±l). 
Therefore, we have 

w<V on9Q+i. (3.19) 

Next, we show that w is a strict subsolution Ciw > in Q^i, if the parameters are chosen 
appropriately. In order to estimate Ciw, we denote 

Ao := - (3.20) 



and notice that 

Wyy = -2x{^x + k) = -2x(/xa; + fc)((l-y^) + y^) = -~2fix{l - y'^){x + Aq) -2ky'^x{-^ 

An 



Then, by a direct calculation and simplification, we obtain 

Ciw = 2^lx{l~y^)h + ky^h, (3.21) 

where 

h - l-2^la{l-y^-Oi + ^~{{b+0:,)-yO^){x + Ao)-y^^^±^0^, 

X X 

h = {l + Oi) + 2iia{l-y^)-2{b+0^)x{^ + l)-2yx{^ + l)Or-,-2y{^ + l)02. 

Aq Aq 

Now we choose r and so that Ciw > holds. Clearly, £iw > if /i, /2 > 0. By (|3.6p . we 
find that, in Q^i, 

h>l-2fia-Cor-{b + N + Cor)Ao, 

^0 

Choose ro to satisfy the smallness assumptions stated above and 

where Co is the constant in (|3.22p . For such a fixed ro, we choose fio to satisfy (|3.17p and 
and Aq to satisfy 



.0 < I, (3.24) 



^^orl<A<^^y (3.25) 

where we have used p.23p to see that 4(70^0 < g^^jj^^y^ in (|3.25p . Then k is defined from 
(13:201) . From (13:221) - (13:251) . 

h, h > 0, 

which implies that 

Ciw > in (3.26) 

whenever r € (0, tq] and /i e (0, fio]. 

By p.l9p . (I3.26p . Remark 1 3. 1[ and the comparison principle f Lemma 13. II) . we have 

^{x, y) > w{x, y) = fix'^{l - y^) - kxy'^ in (5+^. 

In particular, 

^{x,0) > nx^ for X e [0, r]. (3.27) 
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This implies (|3.15p . thus (|3.14p . The proof is completed. □ 
With Proposition 13. 11 we now make the C^'" estimate of ip. 



3.3. C Estimate of xf). If ip satisfies p.2p - p.4l) and (|3.9p . it is expected that is "very 
close" to 1^, which is a solution to (|3.8|) . More precisely, we now prove p.l3p . To achieve 
this, we study the function 

Wix,y):=^~i,{x,y). (3.28) 

By dSIS]), W satisfies 

C2W ■.= ix + aW., + Oi)W,^+02W^y 

Oi - xOa , (3.29) 

W{0,y)^0 on 9Qt^ n {x = 0}, (3.30) 

'^—^x<W,ix,y)<{M+-)x in Qt (3.31) 

a a 

Lemma 3.2. Let a,b, N, R,f, P, and Oi be as in Theorem 13.11 Let fi be the constant 
determined in Proposition 13.11 Then there exist ai G (0, 1) and ri > such that, if 
W e C(Qt^) n C2(g+^) satisfies (13:29)1 - (I33TI1 . then 

W{x,y)<^^x'+- ir^Qtjn/,^ (3.32) 
whenever a G (0, ai] and r G (0,7"!] with fii :— min(2a/i, 1/2). 



Proof. In the proof below, all the constants depend only on the data, i.e., a,b, N, P, R,r, 
Q+^n{x>f/2} V', unle 
By Proposition [3111 



infg+ V'l unless otherwise is stated 



Wix,y)<^—^x^ mQ+.^, (3.33) 



2a ^ro,15i?/16 

where tq depend only on a, 6, R, f, and /3. 
II 
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Fix yo with \yo\ < We now prove that 



Wix,y„)<^^J^x'+^ forxG(0,r). 

By a scaling argument similar to the one in the beginning of proof of Lemma 13. 1[ i.e., 
considering the function ilj{x,y) — ip{x,yo + -^y) in (5^2' conclude that, without loss of 
generality, we can assume that yo = and R = 2. That is, it suffices to prove that 

W(x,0) < i— f^j. ^ ^ X .3 
2ar" 

for some r G {O,ro),a G (0,q;i), under the assumptions that (|3.29p - (|3.3ip hold in Qt^ and 
(133311 holds in Q+ 2- 

For any given r G (0, ro), let 



1 - Ml p 1 - Ml 

, -Dl = 

2a 

z;(a:, y) = ylix2+"(l - j/^^ ^ B^x^y^ (3.36) 



^ Q -"^ ^ Ml n ^ Ml CO Qc:^ 

Air = — , Bi = — , 3.35) 

2a 2a 
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Since (P?^ holds on dQf^ n {x = 0} and (P?^ holds in (5+ 2: then, for all x e (0, r) and 
\y\ < 1, we obtain 

'v{0,y) = = W{0,y), 

v{r, y) = (Air"(l - y^) + B^y^y - > W{r, y), 

^v{x,±l) = Bix^ = > Wix,±l). 

Thus, 

W<v on9Q+i. (3.37) 
We now show that C2V < C2W in From ((X^ . 

£2f - £21^ = £21^ • 

a 

In order to rewrite the right-hand side in a convenient form, we write the term Vyy in the 
expression of C2V as (1 — y'^)vyy + y'^Vyy and use similar expressions for the terms Vxy and 
Vy . Then a direct calculation yields 

jr^y _ ~ ""^^ ^ (2 + - y^)J, + 2Bixy' J2, 



where 

Ji = (l + a)(l + a((2 + a)^ia;"(l-y2) + 2Siy2) -(2 + 04)+ri, 

J2 = 1 + a((2 + - j/^) + 2Biy2) + ^ _ (2 + O4) + T2, 

Ti = (2 + a)!4ixi+" (202a;y(2Si - (2 + a)A^x'^) + 2x2(Si - + O3) + O5 

Oi - a;04 
a 

_ (2 + a)Aix^ 
- 2BI^ 
Thus, in Q^i, 

£2W-/:2VF<0 ifJi, J2<0. (3.38) 

By p.6p and p.35p . we obtain 

iTil, IT2I <Cri-" inQ+i, 

so that, in Q^^, 

< (l + a)(l + ^(l-m)) + (3.39) 

^2 < 1 + ^^(1 -Ml) -2 + Cri-". (3.40) 

Choose ai > 0, depending only on jii, so that, if < a < ai, 

(l + a)(l + ^(l-;,i))-2<-^. (3.41) 

Such a choice of ai > is possible because we have the strict inequality in ()3.41|) when 
a = 0, and the left-hand side is an increasing function of a > (where we have used 
< Ml < 1/2 by reducing ^ if necessary). Now, choosing ri > so that 

ri<min{(^)'^,ro} (3.42) 

is satisfied, we use p.39p - p.4ip to obtain 

Ji, J2 < in Q^-^. 
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Then, by (|3.38p . we obtain 

L2V < C2W in (3.43) 

whenever r S (0, ri] and a E (0, ai]. By (|3.37p . (|3.43p . Remark [3.11 and the standard 
comparison principle (Lemma 13. ip . we obtain 

W <v in Q+i. (3.44) 

In particular, using (|3.35p - (|3.36p with y = 0, we arrive at p.34p . □ 

Using Lemma [3T2l we now generalize the result (|3.32l) for any a € (0, 1). 



Proposition 3.2. Let a, b, N, R, f, (3, and Oi he as in Theorem \3.1[ Then, for any a G (0, 1), 

there exist positive constants r and A which depend only on a, b, N, R, f, f3, and a so that, 
ifW e C(Qt^) n C^iQXn) satisfies (IM) ~ (l33T]) . then 

W{x,y)<Ax^+" m (3.45) 

Proof. As argued before, without loss of generality, we may assume that R = 2 and it 
suffices to show that 

W{x,0) < for X £ [0, r]. (3.46) 

By Lemma [3.21 it suffices to prove (|3.46p for the case a > oi\. Fix any a € (ai, 1) and set 
the following comparison function: 

u(x, y) = ^'^^ x^+"(l - 2/2) + (3 47) 

By Lemma [3.21 

W <u on ag+i for r G (0,ri]. (3.48) 
As in the proof of Lemma 13.21 we write 



Ol - xOi ^ (1-A^i)xl+" ,^ 2w^.O^ , (l-/il)xl+"^ 2f 

= + 2<V^-^. - J/ ) Ji + (2 + "i)^^^^^ 2/ ^2, 

where 



= ^((1 - y^)(2 + «)(^)" + y^(2 + a,)i^^r), 
ii = (l + a)(l+(^)"'i?o) -2 + fi, 



J2 = (1 + ai)(l + i^T'Do) - 2 + f2 



2 



a(l-/xi)):ri+" l^^" " + ""^^"^^ " ^^^^ ^ , 

a + ai)(l - ui)xi+"i V \\ / ' a J 



2ar^i _ , „ , Oi-a;04 

(a + ai)(l - ^1) 
By (|3.6p . we have 

|Ti|,|f2| < Cri-"i 

for some positive constant C depending only on a, 6, iV, /3, ri, and a\. Thus, we have 
max{Ji, J2} < (l + a)(l+(^)"'(2 + a)i^^)-2 + Cri-"i for (x, y) e Q+i. (3.49) 
Choosing r > sufficiently small, depending only on ri, a, and C, we obtain 

Then Lemma |3 . II implies that 



£2^ - li2W ^ C2U - — — < in Q+i. 



W <u in g+i. 
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Thus, ^^M\i holds with 

1-^1 



A 



□ 



Lemma 3.3. Let a, b, N, R, r, (3, and Oi he as in Theorem Vd.ll Then there exist r2 > and 
a2 e (0, 1) such that, if W e C(Qf^) n C^CQ^fl.) satisfies ([3:291) - (f33T|) . we have 

W{x,y)>-^x'+- tnQtj^^, (3.50) 
whenever a G (0, and r £ (0, r2]. 

Proof. By (0301) -(ESD), it can be easily verified that W{x,y) > -^x'^ in Qt^. Now, 
similar to the proof of Leinma l3.21 it suffices to prove that, with the assumption R ~ 2, 

W(x, 0) > -^r-^x^^" for X e (0, r) 

for some r > and a G (0, 02). 

For this, we use the comparison function: 

v{x, y) := -Lx^+"(l - y^) - Kx^y'^, with Lr" = K = i-^. 

Then we follow the same procedure as the proof of Lemma \3l2\ except that C2V > C2W, to 
find that the conditions for the choice of Q;,r > are inequalities (|3.4ip - p.42p with /ii,ri 
replaced by /3, r2 and with an appropriate constant C. □ 

Using Lemma 13.31 we now generalize the result (j3.50p for any a e (0, 1). 

Proposition 3.3. Let a,b, N, R,r, (3, and Oi be as in Theorem 13.11 Then, for any a e 
(0, 1), there exist positive constants r and B depending on a,b, N , R,r, (3, and a so that, if 
W e C(gt^) n C^{QXr) satisfies ^Mf-^Ml, then 

W{x,y)>~Bx'+-^ « Q^3fl/4- (3-51) 

Proof. For fixed a £ (02, 1), we set the comparison function: 



u(x,y) = i ^ - t/) _ L^a;2+a. 2 

Then, using the argument as in the proof of Proposition we can choose r > appropri- 
ately small so that 

C2U- - C2W = £2"- > 

a 

holds for all (.T,?;) e Q+p □ 
With Proposition 13 . 2ft531 we now prove Theorem 3.1. 

3.4. Proof of Theorem 13.11 We divide the proof into four steps. 

Step 1. Let ^ be a solution of p.2p in Qf ^ for r as in Remark [3. 11 and let the assumptions 
of Theorem 13.11 hold. Then ip satisfies p.lOp . Thus it suffices to show that, for any given 
a e (0,1), there exists r > so that 2p e C'^'"(Q^_r,/2) and ^pxx{0,y) = ^, tpxy{0,y) = 
V',;2,(0,y) = 0foraU I2/I < f. 

Let W{x,y) be defined by (|3.28p . Then, in order to prove Theorem 13.11 it suffices to 
show that, for any given a S (0, 1), there exists r > so that 

(ii) D^W{0,y) = for all \y\ < f . 
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Step 2. By definition, W satisfies p.29|) - p.3ip . For any given a G (0,1), there exists 
r > so tliat botli (P^IS)) and (P^Tjl hold in (9+3^/4 by Propositions Fix such r > 0. 

Furthermore, since W satisfies estimate p.3ip . we can introduce a cutoff function into 
the nonhnear term of equation p.29p . i.e., modify the nonlinear term away from the values 
determined by (|3.3ip to make the term bounded in Wx/x. Namely, fix C S C°°(R) satisfying 



_L^<^<Af + - onM; ({s) = s on(-^^, M+-) 
a a a a 

C=0 on M\(-^^, M+-). 
Then, from ([3?29l) and ((33TI1 . it follows that W satisfies 
x{l+aCi^) + ^)Wxx+02Wxy + ib+03)Wyy 



(3.52) 



(3.53) 



Step 3. For z := {x, y) E _r/2' define 



Then 



R.:={M : \s-x\<^,\t~y\<^). (3.54) 

Fix zo = (a;o, j/o) G Qt/2 r/2- Rescale W in Rz^ by defining 

W^^^\S,T)^^W{x^ + '^S,V^ + ^T) for (5,T)egi, (3.56) 

where Qn = i~h, h)^ for h > 0. Then, by ([g^iS)) . (P3T|) . and ([XTO)) . we have 

llM^(-)|ico(gr)<^- (3-57) 
Moreover, since satisfies equation p.29p . W^^°^ satisfies the following equation for {S, T) e 

Qi- 

(1 + f ) (1 + aC(^^^^$^) + 0['°^) W^;f + 0^^° Vi"^") + (6 + 0^^)W^:[,'^^ 

^ + 8 ^ (3.58) 

where 

6['^'\S,T) = . \ g, Oi(^,;/), 0^^(5,T) = -^02(.T,y), 
a;o(l + f) 

4^'')(5,r) = 03(2:,y), d['°\s,T)=Oi{x,y), di^"\s,T) ^ ^05{x,y), 

with a; = a;o(l + f ) and y = ya + ^T. Then, from (I3H)-(I321), we find that, for all 
(5,T) e Qi and zo £ Qr/2.R/2 ^^^^ r <1, 

\di'''\s,T)\<2NV^ for fc = l,...,5, 

\Dd["°\s,T)\ <2NV^ forfc^2, (3.59) 
|L'0^^°^(5,r)| < 27V. 
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Also, denoting the right-hand side of p.58p by F'-''°'' {S,T), we obtain from p.59p that, for 

'r/2,_R/2' 

\F^'°\S,T)\ < Cr^-°', \DF^'°\S,T)\ < Cr^-"", (3.60) 



all {S, T) e Qi and zq e Q+ 



where C depends only on N and a. 
Now, writing equation (|3.58p as 

2 2 

^ Aij{DW^'-''\S,T)Dlw''''>'^ +^Bi{S,T)DiW^'''^ ^ F^'"^ in Qi, (3.61) 

ij = l 2=1 

we get from p. 521) and (|3.58P " (|3.60p that, if r > is sufficiently small, depending only 
on the data, then (|3.6ip is uniformly elliptic with elliptic constants depending only on h 
but independent of zq, and that the coefficients Aij{p,S,T), Bi{S,T), and F^^°^S,T), for 
p e {S,T) e Qi, satisfy 

W^vWcHR^xOT) ^ \\{Bi, ^i/2-a)llci(Qr) - 

where C depends only on the data and is independent of zq. Then, by [101 Theorem Al] 
and (|337| . 

\\W^'"Hc^..,^^^ < C(||W^(-)||^o(Q^) + < C(-i- +rV2-) C, (3.62) 



where C depends only on the data and a in this case. From (|3.62p . 
\DlDlWixo,yo)\ < for aU (x^yo) e Q+/2,,?/2' < * + j < 2. (3.63) 



Step 4. It remains to prove the C"-continuity of D W in Q^^2 /j/2- 
For two distinct points zi = (xi, yi), Z2 = (2:2, 2/2) <= _r/2' consider 

^.^ IW^.x(zi)- W^..(Z2)| 
kl-^2h 

Without loss of generality, assume that xi < X2- There are two cases: 
Case 1. zi e Rz2- Then 



X2 \/ X^, 

xi — X2 + — 5, J/1 = 2/2 H —T for some (5, T) e Qi. 



By (EH, 

\w'^sf{S,T)-W'^,f{OM ^ 

(52 + T2)a/2 - ' 

which is 



< c. 



< C. 



((xi -a;2)2+a;2(yi-y2)2)"/2 

Since X2 G (O,?") and r < 1, the last estimate implies 

\Wxxixi,yi) - Wxx{x2,y2)\ 

{{xi - X2f + {yi - y2fr/^ 

Case 2. zi ^ R^^. Then, either |xi - 2;2| > ^ or \yi - 2/2! > Since < a;2 < r < 1, 

we find 

i--..r>(f)°. 

Thus, using (|3.63p and xi < X2, we obtain 

\Wxxizi)-Wxx{z2)\ ^ \Wxxizi)\ + \Wxx{z2)\ ^ ^ xf + ^ 

Ui-Z2|" " U1-Z2I" - a;? - ■ 
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Therefore, A < 2C in both cases, where C depends on a, r, and the data. Since zi ^ Z2 
are arbitrary points of Q^/j ji/2: obtain 




(3.64) 
cs, we 






where C depends on a,r, and the data, but is independent of J > and zq. 



Thus, W G C^'°'[Qt RI2) ^^^^ II ^1^2 ar7)+ depending only on the data because r > 

depends on the data. Moreover, p.63p imphes D^W{0,y) = for any \y\ < R/2. This 
concludes the proof of Theorem 13.11 □ 

4. Optimal Regularity of Solutions to Regular Shock Reflection across 



As we indicated in Section 2, the global solution cp constructed in [10] is at least C^'^ 
near the sonic circle P1P4. On the other hand, the behavior of solutions to regular shock 
reflection has not been understood completely; so it is essential to understand first the 
regularity of regular reflection solutions. In this section, we prove that C^'^ is in fact the 
optimal regularity of any solution ip across the sonic circle P1P4 in the class of standard 
regular reflection solutions. Our main results include the following three ingredients: 

(i) There is no a regular reflection solution that is across the sonic circle; 

(ii) For the solutions constructed in [TU] or, more generally, for any regular reflection 
solution satisfying properties ([u]) and ([iv])~(Ivil) at the end of Section 2, ip is C^'" in the 
subsonic region fl up to the sonic circle P1P4, excluding the endpoint Pi, but D'^ip has a 
jump across P1P4; 

(iii) In addition, D^ip does not have a limit at Pi from fl. 

In order to state these results, we flrst deflne the class of regular reflection solutions. As 
proved in [10| . when the wedge angle 0„ is large, such a regular reflection conflguration exists; 
and in 111], we extend this to other wedge-angles for which regular reflection conflguration 
exists. 

Now we deflne the class of regular reflection solutions. 

Definition 4.1. Let 7 > 1, pi > po > 0; ^.i^d, ui > be constants, and let be defined by 
(|2.15[) . Let the incident shock 5 = = ^0} hit the wedge at the point Pq = (^Oi tan^^,), 
and let state (0) and state (1) ahead of and behind S be given by ()2.13p and (|2.14p . respec- 
tively. The function ip € C'^'^{A) n C^{fl) is a regular reflection solution if ip is a solution 
to Problem 2 satisfying p. 2111 and such that 

(a) there exists state (2) of the form (|2.18p with U2 > 0, satisfying the entropy condition 
P2 > Pi and the Rankine-Hugoniot condition {piDipi — p2Dip2) ■ v — along the line 
Si :— {ipi — ip2} which contains the points Pq and Pi, such that Pi £ A is on the sonic 
circle of state (2), and state (2) is supersonic along PqPi; 

(b) equation p.7p is elliptic in fl; 

(c) ip > ip2 on the part P1P2 = ^shock of the reflected shock. 

Remark 4.1. If state (2) exists and supersonic, then the line Si = {ipi — ip2} necessarily 
intersects the sonic circle of state (2); see the argument in [10] starting from (3.5) there. 
Thus the only assumption regarding the point Pi is that Si intersects the sonic circle within 



THE Sonic Circle 



A. 



REGULARITY OF SOLUTIONS TO REGULAR SHOCK REFLECTION 



19 



Remark 4.2. The global solution constructed in |10| is a regular reflection solution, which 
is a part of the assertions at the end of Section 2. 

Remark 4.3. There may exist a global regular reflection configuration when state (2) is 
subsonic which is a very narrow regime [14U36j . Such a case does not involve the difficulty 
of elliptic degeneracy which we are facing in our case. 

Remark 4.4. Since tp — ipi on Tshock by (j2.2ip . Condition (c) in Definition \A.l\ is equivalent 
to 

Tshock C {(p2 < Vl}, 

that is, Tshock is below Si. 

Furthermore, we have 
Lemma 4.1. For any regular reflection solution Lp in the sense of Definition 14. ll 

if > Lp2 in f2. (4-1) 
Proof. By ^H^-^^ and ([^1^ . i) := Lp - Lp^ satisfies 



(c^ - (-0? - C + W2)^)'0?5 + (c^ - (V'r, ~V + U2 tan6li„)^)V',,,, 

— 2{7p^ — C + U2)(V-'») — ?? + 1*2 ta,n9iu)ip^r] = in Q, 



(4.2) 



where c^{D'ip,->p,^,ri) = + (7 - - U2)ip^ + iv - U2 tan6'u,)'0>7 - ^\D^? - i' 

Since the equation (|4.2[) is elhptic inside and (p is smooth inside fi, it follows that (|4.2p 
is uniformly elliptic in any compact subset of Q. Furthermore, we have 

■0 = on Tsonic, 

Dip ■ (- sin6'^, cos6'u,) on T^edgc, 

i/",, = — W2 tan^^ < on dfl n [i] = 0} , 

■0 > on Tshock (by Definition 4.1 (c)). 

Then the strong maximum principle implies 

^ > in 

which is (j4.1|) . This completes the proof. □ 



Now we first show that any regular reflection solution in our case cannot be across 
the sonic circle T sonic ■= PiPa- 

Theorem 4.1. Let ip be a regular reflection solution in the sense of Definition 4.1. Then 
p cannot be across the sonic circle Ts onic ■ 

Proof. On the contrary, assume that ip is across T sonic- Then i/j = ip — p2 also 
across Tgonic, where (p2 is given by (|2.18p . Moreover, since tp = in P0P1P4 by (|2.21|) . we 
have DV(C,^) = at all (C,?/) € Ts onic • 

Now substituting ip = ip + p2 into equation (|2.7|) and writing the resulting equation in 
the (a:, ?/)-coordinates (|2.20p in the domain ilgj, defined by ()2.24p . we find by an explicit 
calculation that ip{x,y) satisfies equation p.2p in with a = 7 + 1 and b = and with 



20 



MYOUNGJEAN BAE, GUI-QIANG CHEN, AND MIKHAIL FELDMAN 



(4.3) 



O-j = 0-i{x,y,il:,tp^,ipy),j = 1, ... ,5, given by 

C2 2c2 C2 2(c2 - a;)"^ 

2 

02{Dtlj,^,x) = ;^ 75" ("02; + C2 - 

C2(C2 - X)^ 

03{Dij,i;,x) = — ^ -(x{2c2-x)-{j-l){i: + {c2-x)i;, + lijl) 

2{c2~xY^y)' 

04{Dij, ^, x) = (x - ^ {i/j + (C2 - x)yj^ + + )) , 

C2 — X\ C2 2 2[C2 — X)'" / 

05{D^,^P,X) = , \ (V'x+2c2-2x)V'y. 

C2\C2 X) 

Let (0, j/o) be a point in the relative interior of Tgonic- Then (0, yp) + Q^r ^eo if 
r, i? > are sufficiently small. By shifting the coordinates {x,y) {x,y — t/o)i we can 
assume (0, yo) = (0, 0) and C i^eo- Note that the shifting coordinates in the y-direction 
does not change the expressions in (|4.3p . 

Since ip G C'^ifleg U T sonic) with D'^ip = on T sonic, reducing r if necessary, we get 
\Dtp\ < 6x in Q^fi, where J > is so small that p.9p holds in Q^^, with (3 = AI = 1, and 
that the terms Oi defined by (|4.3p satisfy (|3.6p - p.7p with M — 1. Also, from Definition 
14.11 we obtain that ?/; = 1^9 — 1^2 > in Q^j^- Now we can apply Proposition 13 . II to conclude 

'4^{x,y) > fix"^ on Q+^5^^^g 

for some /i,r > 0. This contradicts the fact that D^il;{0,y) = for all y e (— i?, i?), that is, 
77) = at any (^, 77) G T^on.c. □ 

In the following theorem, we study more detailed regularity of "0 near the sonic circle 
in the case of C^'^ regular reflection solutions. Note that this class of solutions especially 
includes the solutions constructed in |10j . 

Theorem 4.2. Let ip he a regular reflection solution in the sense 0/ Definition 14.11 and 
satisfy the properties: 

(a) ip is C^'^ across the part T sonic of the sonic circle, i.e., there exists Eq > such that 
ip e Ci'H^WWf^n {c2 - eo < r < C2 + eo}); 

(b) there exists (5o > so that, in the coordinates (|2.20|) . 

2 (5 

\dx{ip - ip2){x,y)\< — -^x in n,g] (4.4) 
7 + 1 

(c) there exist uj > and a function y ~ f{x) such that, in the coordinates (j2.20p . 

n,,={{x,y):xe{0,eo),0<y<f{x)}, 

(4.5) 

^shock n {0 < X < £0} = {{x, y) : xe (0, eo), y = f{x)}, 



J r 

ll/llci.i([o,so]) < 00 , —>w>0 /or 0<a;<eo- (4.6) 

Then we have 



REGULARITY OF SOLUTIONS TO REGULAR SHOCK REFLECTION 



21 



(i) Lp is C up to T sonic away from the point Pi for any a G (0, 1). That is, for any 
a € (0, 1) and any given (^Oi ^o) G ^sonic \ {^'i}; there exists K < oo depending only 
onpo,pi,"f, So, a,\\ip\\ci.i{n,g), and d dist{{^o,T]o), Tshock) so that 

ll'^ll2,Q;Sd/2({o,r/o)nn,o/2 - 

(ii) For any (^o,'7o) G T^soriic \ {Pi}, 

1 



lira {Drr'^ — Drr^2) 



(S>»7)->(«0.'I0) 7+1 

(iii) D'^ip has a jump across Tsonic-' For any (Co,'?o) G ^sonic \ {Pi}, 

lim Drrf — lim Drr^ = ; 

(S.'7)^(«0.')0) (S.v)^{io,Vo) 7+1 

(5,»j)efi (5,T,)eA\s-! 

(iv) The limit lim(5,,,)^f.i D^ip does not exist. 

Proof. The proof consists of seven steps. 
Step 1. Let 

By (|2?2T|) and (144]) . we have 

V'(0,y) = 7/;,(0,y) = Vy(0,2/) = for all (0, j/) € T.^^e, (4.7) 
and thus, using also (14.511 - (14.611 . we find that 

|V(a;,y)| <Cx2, \D,.yijix,y)\ < Cx for all (x, y) e f^^o, (4.8) 

where C depends only on ||V'llci.i(n;;;) ^^^^ ll/llci([o, £o])- 

Recall that, in the {x, j/)-coordinates p.20p in the domain ^1^^ defined by (|4.5p . ip{x, y) sat- 
isfies equation p.2p with Oi = Oi{x,y,ip,ilJx,il'y) given by (|4.3p . Then it follows from (|4.3p 
and (|4.8p that (|3.6p - (|3.7p hold with TV depending only oneo, ll^llci i(ni^)' ^^'^ ll/lici([o, eq])- 

Step 2. Now, using (14. 4p and reducing eq if necessary, we conclude that p.2p is uniformly 
elliptic on Ci {x > 6} for any 5 e (0, £o)- Moreover, by (c), equation ()3.2p with ()4.3p . 
considered as a linear elliptic equation, has coefficients. Furthermore, since the boundary 
conditions (|2.16p hold for ip and ip2, especially on Tyjedge = {u = 0}, it follows that, in the 
(x, ?;)-coordinates, we have 

'ipy{x,0) = for all X G (0,eo)- (4.9) 

Then, by the standard regularity theory for the oblique derivative problem for linear, uni- 
formly elliptic equations, ip is in fJe^ up to dileg D {0 < x < Eq, y = 0}. From this and 
(c), we have 

^ e c^'H"^,) n c\n,, u rt"]^,), (4.io) 

where r^^^"]^^ T^edge n {0 < a; < £o} = {(a;, 0) : < a; < Eq}- 
Refiect l^eo with respect to the y-axis, i.e., using (|4.5p . define 

n,, {{x, y):x& (0, Eq), ~fix) < y < fix)}. (4.11) 

Extend ip{x,y) from fJ^g to fJ^^ by the even reflection, i.e., defining ^(x,—y) — ili{x,y) for 
{x,y) G rigQ. Using (|4.9p - (|4.10p . we conclude that the extended function il;{x,y) satisfies 

7/.eCi^i(a^)nc2(A,J. (4.12) 

Now we use the explicit expressions p.2ll and (|4.3p to find that, if ip{x,y) satisfies equa- 
tion p.2p with (|4.3p in Og^, then the function ip{x,y) :— ■(/'(a;, — j/) also satisfies p.2p with 
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Ok{Dip,ip, x) defined by (|4.3p in ^1^0- Thus, in the extended domain Ct^^, the extended 
tp{x, y) satisfies p.2[) with Oi, . . . , O5 defined by the expressions (|4.3p in Cl^^ . 

Moreover, by (|2.2ip . it foUows that -0 = on T sonic- Thus, in the (x, y)-coordinates, for 
the extended "0, we obtain 

0(0, v) = for aU y e (-/(O), /(O)). (4.13) 

Also, using > (/32 in r2, 

V'(0,2/)>0 in (l,,. (4.14) 

Step 3. Let P = (^*,77*) e Fsonic \ {^i}- Then, in the (x, y)-coordinates, P = (0,y*) 
with j/» e [0,/(0)). Then, by (|4.6p and (|4.1ip . there exist r,R > 0, depending only on 
eo,C2 C2(Po,Pi,"i,^«)), and d = dist((^* , 77*), r^hocfe), such that 

(0,2;*) + Q+fl, c fieo- 

Then, in Q^^, the function ip{x^y) :— ipix,y — y*) satisfies all the conditions of Theorem 
13.11 Thus, applying Theorem 13. II and expressing the results in terms of -0, we obtain that, 
for aU y, E [0,/(0)), 

<.Ji?o,„,, '/'-(^'y) = r^' (.JiTo.«.) ^-y^'^'y^ = iJ^^a.y,, ^vvi^^y) = »• (4-i5) 



Since 0rr = '(pxx by (|2.20p . this implies assertions of Theorem 14.21 

Now assertion (jml of Theorem l4.2l follows from (jHI since, by ()2.2ip . (/j = 1^92 in i3e(C*: ??*)\^^ 
for small £ > and if 2 is a C°°-smooth function in R-^. 

Step 4. It remains to show assertion ([iv| of Theorem l4.2l We prove this by contradiction. 
Assume that assertion (jlvl is false, i.e., there exists a limit of D^0 at Pi from fi. Then our 
strategy is to choose two different sequences of points converging to Pi and show that the 
limits of 4'xx along the two sequences are different, which reaches to a contradiction. We 
note that, in the (a;, y)-coordinates, the point Pi = (0, /(O)). 

Step 5. A sequence close to Tsonic- Let {y'"m}m=i be a sequence such that j/m' e 
(0, /(O)) and hm,„^oo = /(O). By there exists x'i^ G (0, i) such that 

|V...(xW,ya)) _ -1-1 + |0.,(a;«,yW)| + < ^ 



for each to = 1,2,3,.... Moreover, using (|4.6p . we have 

< /(O) < /(x,W). 

Thus, using (|4.5p . we have 
(a;W,yW)el], hm (x«,£)) = (0,/(0)), 



m — >-oo 



1 (4-16) 

lim 0^:r(a;^\yii') = — — , lim 0a.y(a;^^ yi^^) = lim 0j;y(a;(;[), y£^) = 0. 

m — *C30 ^ + i — *oo rn — »oo 

Step 6. The Rankine-Hugoniot conditions on T shock- In order to construct another 
sequence, we first combine the Rankine-Hugoniot conditions on T shock into a condition of 
the following form: 

Lemma 4.2. There exists e € (0,eo) such that -0 satisfies 

bi{x, y)il:x + hix, y)ipy + bs^x, y)ip = on Tshock n {0 < a:; < e}, (4-17) 



where bk G C {Tshock H {0 < a; < e}) and further satisfies 

bi{x,y)>\ |&2(a;,y)| < |63(a;,y)|<^ on T shock ^ {'i < x < e} (4.18) 
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for some constant A > 0. 

Proof. To prove this, we first work in the ry)-coordinates. Since 

ip ^ (fii, pDip ■ V = piDipi ■ V onTshock, 
then V is parallel to Dipi ~ Dlp so that 

{piD^i - pDip) ■ {Dipi ^Dip)=0 on Tshock- (4.19) 
Since both cp and Lp2 satisfy (|2.8p - (|2.9p and 'Z' — V2, we have 

-i 1 

p = p{D^,^b,t V) = {pV + (7 - - "2)V'« + iv- V2)4'n - - ^)) ' ' ' 

(4.20) 

Then, writing if = (p2 + ip and using (|2.13p - (|2.14p . we rewrite (|4.19p as 

E{tP^,tP^, V, ^,V)=0 on Tshock, (4.21) 
where, for (pi,P2,P3, ??) € 1^^- 

S(pi,P2,P3,e,^) (4.22) 

= Pl(("l - 0(W1 - "2 -Pl) + f?(w2 +P2)) 

-p{Pl,P2,P3,i,V)(^iu2 -^+Pl){ui - U2 -Pl) - (t'2 - ?7+P2)(v2 +P2)), 

P(Pl,P2,P3,e,^) = (pr' + (7- l)((e - W2)P1 + V2)P2 ~ \\P\^-P3)) '"'(4.23) 

with V2 '■= U2 tan But- 

Since both points Pq and Pi lie on Si = {ipi = ^2}, we have 

{ui - U2){£.i - Co) - U2('7i - ^70) = 0, 

where {^i,rii) are the coordinates of Pi. Now, using the condition ip = ip i on Tshock, i.e., 
ip + ip2 ^ ipi on Tshock, we have 

(Ml -M2)(C~Cl) - V'(C,^?) „ -p 

?/ = h r^i on Ts/iocfc. (4.24) 

V2 

From (|4.2ip and (|4.24p . we conclude 

F(i/;5,V^,V,C) =0 onr,,,o,fc, (4.25) 

where 

i?r £\ T?( c (wi - U2)(C - Ci) -P3 , X . . „„x 

F{pi,P2,P3,0 = E{pi,p2,P3,^, h?7i)- (4.26) 

V2 

Now, from (|3221)-(|1!2S1), we obtain that, for any ^ G R, 
FiO, 0, 0, = EiO, 0, 0, t ^"^'"^^^^'^^^ + m) 

1^2 

= Pl(^{ui - ^l){ui - U2) + V2lJ^ - P2[{U2 - 6)("1 ~ U2) - V2{V2 - Tji)) 

= [piDipi{(,i,r]i) - p2Dip2{£,i,-qi)) ■ {ui - U2, -V2) 
= 0, 

(4.27) 

where the last expression is zero since it represents the right-hand side of the Rankine- 
Hugoniot condition (|2.12p at the point Pi of the shock Si = {pi — 1P2} separating state (2) 
from state (1). 
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Now we write condition (I4.25P in the (x, y)-coordinates on Tghock n {0 < x < Eq}- By 
(I2HD-(I22DD and (112S1), we have 

^{^^,^py,'il;,x,y) = on r^^ocfe n {0 < a; < eo}, (4.28) 

where 

'^{Pi,P2,P3,x,y) ^F(-pi cos{y + 9^) sin(y + 6'„), 

- pi sm{y + e.u,) -\ cos(y + 6*^,), ps, U2 + (c2 - x) cos(y + 9^)). 

0-2 — X 

(4.29) 

From (|4:27l) and (jiTM)) . we find 

^'(0,0,0,x,j/) = F{0,Q,0,U2 + {c2-x)cos{y+d^)) = on r,hockn{Q <x< Eq}. (4.30) 

By its exphcit definition ([iT^ - fi?^ . and the function *(pi,p2,P3, a:, ?/) 

is C°° on the set {\{pi,p2,P3,x)\ < S}, where S > depends only on U2,V2, P2, ^o,Vo^ i-e., 
on the data. Using (14. 8p and choosing e > small, we obtain 

\x\ + \'iP{x,y)\ + \DtP{x,y)\ <S for aU {x,y) G H;. 

Thus, from it follows that -0 satisfies (|IT7)) on T shock n {0 < a; < e}, where 

bk{x,y)=J 'i'p^(^til;^{x,y),tiljy{x,y),tij{x,y),x,y^ dt forfeit, 2, 3. (4.31) 
Thus, we have 

1 



bk e C{Tshock n {0 < a; < e}), |6fc| < - on Tshock n {0 < a; < e}, for fc = 1, 2, 3, 
for some A > 0. 

It remains to show that 6i > A for some A > 0. For that, since bi is defined by (|4.3ip . we 
first show that ^Pp^ (0, 0, 0, 0, yi) > 0, where {xi,yi) = (0, /(O)) are the coordinates of Pi. 

In the calculation, we will use that, since (0, yi) are the (a;, y)-coordinates of Pi = (^i, 771), 
then, by flJ ^ - fi:^ . 

^1 =U2 + C2 cos(2/i + 9ni), Vi = ^2 + C2 sin(2/i + 9w), 

which implies 

(6 - U2f + {m - W2)' = cl 

Also, c| = P2~^- Then, by explicit calculation, we obtain 

(0, 0, 0, 0, yi) ^^{{ui - ei)(ei - "2) - m im - V2)) 

P2 (^-^^^ 

{{U2 - Cl)(6 - W2) + {V2 - '7l)('7l - '"2))- 

C2 

Now, working in the (^, 77)-coordinates on the right-hand side and noting that Dipi (^i , 771 ) = 
(ui -£,i,-rji) and Dip2{£,i,'ni) = {u2 - ^1 , V2 ~ rji) , we rewrite (|4.32|) as 

^p, (0, 0, 0, 0, yi) = --^ (pii?^i(ei, r?i) - p2£'</'2(ei, ^1)) ' i?</?2(6, ?7i), 

where D = (9^,i9^). Since the point Pi lies on the shock 5*1 = {tpi = 1P2} separating state 
(2) from state (1), then, denoting by tq the unit vector along the line Si, we have 

D'Pi{^i,m) ■ Tq = Dip2{^i,m) ■ To- 
Now, using the Rankinc-Hugoniot condition (|2.12p at the point Pi for ipi and (^2, we obtain 

piD(pi{£,i,m) - P2D(p2i£.i,m) = {pi - P2){Dif2{^i,m) ■ to)to, 
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and thus 

where p2 > pi by the assumption of our theorem. 

Thus it remains to prove that Dip2{^i,Vi) ' '^o 7^ 0- Note that |D(p2(Ci, '7i)| — ^2 = 
p^2 since (Cij^yi) is on the sonic circle. Thus, on the contrary, if Dip2{£,i,'ni) ■'^0 = 0, 

then, using also ryi) • tq = iD(^2(Ci,f?i) • tq, we can write the Rankine-Hugoniot 

condition (|2.12p at {£,i,rji) in the form: 

p,\D^,{^,,^,)\^ P2P^r'^'' = p^r'~"'- (4.33) 

Since both Lpi and (^2 satisfy p.7p and since <y3i(^i,?7i) = <(52(Ci,f?i) and \Dip2{^i,rii)\ — C2, 
we have 

7-1 I ~ I n m2 7-1 ,7^1 7-1 

Pi + =P2 +—3—^2 • 

Combining this with (j4.33p . we obtain 

2 

= 1. (4.34) 



2 

7 + 

Consider the function 







■1' 






7 + 







2 .7-1^1^.-2 



g(s) = -—s-^-' + ^s-^ on (0, 00). 
7+1 7+1 

Since 7 > 1, we have 

5'(s)<0 on (0,1); g'(s) > on(l,oo); .9(1) = 1. 

Thus, g{s) = 1 only for s = 1. Therefore, (j4.34p implies pi = p2, which contradicts 
the assumption pi < p2 of our theorem. This implies that D(/52(^i, f?i) • t'o 7^ 0, thus 
*p,(0,0,0,0,yi) >0. 

Choose A := ^^fp^ (0, 0, 0, 0, j/i). Then A > 0. Since the function '^{pi,p2,P3,x,y) is C°° 
on the set {\{pi,p2,P3, x)\ < S} and since ^/i G C^'^(n^) with V'(0,0) = ^'2^(0,0) = i/'iy(0,0) = 
by (|4.7p . we find that, for small e > 0, 



^p^{t%j;x{x,y),t'4!y{x,y),til:{x,y),x,yj > A for aU (x,?/) G Tshock(^{0 < a; < e}, t <E [0,1]. 
Thus, from (|4.3ip . we find bi > A. Lemma l4?2] is proved. □ 

Step 7. A sequence close to T shock- Now we construct the sequence close to Tghock- 
Recall that we have assumed that assertion (jrv)) is false, i.e., D^^/; has a limit at Pi from 17. 
Then (jiTTB)) implies 

lim ^xyix,y) ^ lim ipyyl^x^y) = 0, (4.35) 

(x,y)^(0, /(O)) (rE,y)^(0, /(O)) 

(a;,H)en (x,y)efi 

where (0, /(O)) are the coordinates of Pi in the (a;, y)-plane. Note that, from (|4.7p . 

ijy{xj{x))^ ij^y{sJ{0))dS + 'ipyy{x,t)dt, 

Jo Jf{0) 
and, from (|2.24p . all points in the paths of integration are within fl. Furthermore, by (|2.25p . 
< f{x) - /(O) < Cx with C independent of a; G (0,£o)- Now, ((135)1 impfies 

lim M^lM = 0. (4.36) 

x—>0+ X 
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Also, by Lemma 

\M^J{x))\ = \b^y±M\ < + 1^1) on (0,8), 

bi 

where e > is from Lemma [4.21 Then, using (|4.36p and \ilj(x,y)\ < Cx^ by (|4.8p . we have 

Um ^iMM ^ 0. (4.37) 

2:^0+ X 

Let 

T{x) := i/j-^Xx, fix) - —x) 

for some constant uj > 0. Then .F(x) is well-defined and differentiable for < a; < Eq so 
that 

:F{x) = tp^ix, f{x) - ^x) 

= /(a;)) + / -^^^(^^ f(x) - Y^x) dt (4.38) 

a; /'"'^ tuj 
= tpx{x, f{x)) - —X J Tpxy{x, fix) - —a;) dt. 



Now (|4?35)) and (|4^ imply 



lim = 0. (4.39) 

a:— »0+ X 



By (OU)) and since / e Ci^i([0, Eq]), we have 

J^e C([0,£])nCi((0,e)). (4.40) 

Then (j4.39p and the mean- value theorem imply that there exists a sequence {a;^^^} with 
x'i^^ e (0, e) and 

limxf^=0 and lim J^'(x[?^') = 0. (4.41) 

k — *oo k — -'Oo 

By definition of Tix), 

ij^^{x, g{x)) = T'{x) - g'ix)i^^yix, g{x)) (4.42) 

where g{x) := f{x) - j^x. 

On the other hand, |/'(a;)| is bounded. Then, using (li?^ and (|0T|) - (|0^ yields 

lim V..(4'^5(4'^)) = 1™ ^'(4'^) =0. 

k — *oo k — >oo 

Note that lim2;^o+ = /(O)- Thus, denoting = g{x\, ), we conclude 

(4^\yf ) e hm (xf ,yf ) = (0,/(0)), hm ^..(4'\2;f)) = 0. 

/C — *C30 /C — *C30 

Combining this with (|4.16[) . we conclude that does not have a limit at Pi from fJ, which 
implies assertion ((iv|) . This completes the proof of Theorem 14.21 □ 

Remark 4.5. For the isothermal case, 7=1, there exists a global regular reflection solution 
in the sense of Definition WA\ when 0^ G (0, §) is close to ^. Moreover, the solution has 
the same properties stated in Theorem 14.21 with 7=1. This can be verified by the limiting 
properties of the solutions for the isentropic case when 7 — > 1+. This is because, when 
7^1+, 

iip) ^ Inp, pip) ^ p, c^ip)^! in 
pilD^l^f) ^ poe^^^+^l'^^l') in (EH) 
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and 

in which case the arguments for establishing Theorem 14.21 is even simpler. 
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